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Abstract. Dissipative processes in physics are usually associated with non-unitary actions. 
However, the important resource of entanglement is not invariant under general unitary 
transformations, and is thus susceptible to unitary "dissipation". In this note we discuss both 
unitary and non-unitary dissipative processes, showing that the former is ultimately of value, 
since reversible, and enables the production of entanglement; while even in the presence of 
the latter, more conventional non-unitary and non-reversible, process there exist nonetheless 
invariant entangled states. 



1. Introduction to Bipartite Entanglement 

1.1. Definition and Measure (Concurrence) 

Bipartite entanglement involves the direct product space V of two (complex) vector spaces, 
V\ and V 2 , of dimension m and n respectively. If V\ has basis {e;, i = 1 . . . m} and V 2 has 
basis {fj, j — 1, . . . n} then {e^ eg) fj, i — 1 . . . m, j — 1 . . . n} is a basis for V — V\ <g) V 2 . 

In this note we specialize to the specific and more familiar case of two two-qubit spaces, 
m = n = 2, with the same standard bases {e,(= f),i = 1,2} for V\ and V 2 . We shall 
also use the matrix forms {ei, e 2 } = {[1, 0], [0, 1]} as well as the ket notation ei = |0) and 
e 1 (g> e 1 = |0, 0) etc. Note that in the context of quantum mechanics we refer to vectors as 
pure states. 

Definition 1 (Entangled pure state) Every vector in V x <g> V 2 is a sum of products; but not 
every vector is a product. If it is a product, then it is said to be non-entangled or separable. 

It is a straightforward matter to determine whether a vector v E V is entangled or not. 
IfveV is non-entangled, i.e. separable, then 

2 

= {x x e x + x 2 e 2 ) g> (t/iei + y 2 e 2 ) 



Vi 0,i = i,2} (i) 



from which we deduce that the matrix c of coefficients c^- has determinant zero, det c = 0. 
\ For typographical simplicity we write all our (column) vectors as row vectors. 
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Example 1 (Separable pure state) Consider the bipartite pure state 

(l/v / 50)(|0, 0> + 3|0, 1) + 2|1,0) + 6|1, 1)). 
The matrix c of coefficients is given by 

c= 1/V50 

for which det c = and so the state is separable. 



1 

2 



(2) 



(3) 



Example 2 (Bell state) An example of a maximally entangled two-qubit state is given by the 
Bell state ^=(|0,0) + |1, 1)) for which 



1 

V2 



1 
1 



(4) 



and so det c = 1/2. 



This simple criterion for pure state separability in fact gives a measure of entanglement for 
pure states. To obtain this measure, we normalize so that a Bell state, such as that of Example 
121 has maximal measure of entanglement equal to 1, and we arrive at the definition of the 
concurrence applicable to pure states: 

Definition 2 (Pure state Concurrence) A measure of entanglement for pure bipartite states 
(belonging to two two-qubit spaces) is given by the concurrence C = 2| det c\. 

The concurrence C is essentially equivalent to the measure S called Entanglement of 
Formation, based on the Von Neumann entropy of the partial trace[l J§|. 

The entanglement measure C may be extended to general, or mixed, states (density 
matrices). We describe this extension in Section |4~3| 



1.2. General states 

We are initially concerned in this note with pure states; i.e. represented by vectors v in V. 
However, we can equally represent our vector v by the matrix p[v) = vv^. Of course the 
overall irrelevant phase information is lost in this form. It is easily verified that p is a hermitian 
matrix of rank one, that is, all sub-matrices of order 2 or more have determinant zero. And it 
has a sole non-zero eigenvalue which is equal to 1. It has trace equal to one, assuming that v 
is normalized. A hermitian matrix all of whose eigenvalues are > is called a positive matrix 
(more accurately, semi-positive). We may extend this description of the matrix associated 
with a pure state to give the following definition of a state in general (mixed state or density 
matrix): 

Definition 3 (State) A state p (acting on a space V) is a positive matrix of trace 1. 
Equivalently, 

Definition 4 (State as convex sum) A state p is a convex sum Y.i \p% {\ > J2i A* = 1) 

of pure states p^. 

We simply note here the definition of separability for (general) states: 

Definition 5 (Separable state) The state p acting on V\ (g) V 2 is said to be separable if is given 
by a convex sum J2i \p\ ® Pi {\ > J2i K — 1) where pf acts on V a . 



§ Writing / = -xlog 2 x - (1 - x)log 2 (l - x) then £ = /((l + Vl -C 2 )/2). 
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When p = p 1 eg) p 2 it is said to be simply separable. The above definition extends immediately 
to multipartite states. 

If we have a measure of entanglement £ for pure states (such as that given in Definition^ 
we may extend it to general states by 

Definition 6 (Entanglement of general state) The entanglement £(p) of a mixed bipartite 
state p acting on V\ <8> V 2 is given by 8{p) = min{J2i P^i^^lp = J2i PiPi^i)} where the ipi 
are pure states in V\ <S> V 2 . 

1.3. Unitary and Local Unitary Transformations 

Since every (normed) vector v E V can be transformed to the (non-entangled) state 1 0, 0) 
by a unitary transformation, it is clear that entanglement is not invariant under unitary 
transformations. However, under a local unitary transformation, defined by U = U\ ® Ui, 
one can see that the concurrence, for example, is invariant: 

Theorem 1 The concurrence C is invariant under local unitary transformations. 

Let v = J2i,j=i...2 a ij e i ® E V = V\ ® V 2 , and the unitary matrix U = U\ <g> U 2 be a local 

unitary matrix; then 

Uv = dijUxti <g> U 2 ej 

where ct r = Q>ij(Ux)ik(U2)jr so that c = U\aU 2 whence 

| detc| = | det(Z7iaZ7 2 )| 

= | det Ui | | det a| | det U 2 \ 

= | det a\ since | det Ui\ = 1. 

We may see rather more immediately from Definition [5] that the property of being 
separable is invariant under local unitary transformations; and this extends to the multipartite 
case. However, an extension of Theorem Q] to multipartite systems, namely that such local 
transformations preserve the measure of entanglement, would depend on a definition of 
measure (or measures) of entanglement for such systems, which is currently unavailable. For 
general multipartite states, local unitary equivalence does not preserve all the relevant (state 
and substate) entanglement properties 

2. Unitary Dissipation 

Although the notion of dissipation is more usually associated with a non-unitary process, 
from the preceding we see that entanglement is subject to unitary dissipation, since 
unitary evolution associated with a (hermitian) hamiltonian does not necessarily preserve 
entanglement. Of course, the good news is the other side of this coin; that is, entanglement 
may be produced by the evolution induced by a quantum control hamiltonian. Quantum 
control applied to multipartite systems has been well treated, see for example [3 |. We choose 
a simple example to illustrate the bipartite case. 
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Example 3 (Entanglement production and decay) Consider the unitary evolution U(t) 
induced by the hamiltonian H given by 



H 











y 





x 2 














x 3 





y 








Xi 



(5) 



This is essentially a free hamiltonian with the addition of an off-diagonal time-independent 
control term Jjjl Note that without this latter term H would not change the entanglement since 
H would then be a local transformation. 

We act by the unitary evolution matrix U (t) = exp(itH) induced by Eq.^ on the base 
vector v = [1, 0, 0, 0]. (Note that for such calculations it is important to choose a fixed basis 
- here the standard basis.) 



v(t) = exp (itH)vo 



e itXl [cos(ty) ,0,0,isin(tj/)]. 



(6) 
(7) 



Apart from an overall phase factor, only the control term y plays a role in the entanglement 
production. 

Using the measure of entanglement for pure states given in Definition^ the concurrence 
forv(t) is given by \ sm(2ty) \ (See Figure 1 ). 



2.1. Unitary dissipation of Entanglement 

Example 4 (Unitary dissipation of entanglement) Referring to the previous example and 
Figure 1, we see immediately that at t = 7r/4 (in units ofl/y where y is the control frequency) 
we have the maximally entangled (Bell) state ^[1,0, 0,z] (up to an overall phase factor). The 
unitary action U(t) destroys the entanglement, completely att = ir/2. 



3. Environmental Dissipation 

For the usual description of dissipative processes we must use the standard definition of a 
general quantum state p given by Definition [3j Thus p is an iV x N positive matrix (and for 
our two-qubit examples, iV = 4). For a non-dissipative process, the basic equation which 
determines the evolution of a hamiltonian quantum system may be written in the form of a 
differential equation for the quantum state, the Liouville-Von Neumann equationflU (choosing 
units in which Ti — 1): 

ij t p(t) = [H,p(t)]=H P (t)-p(t)H (8) 

where H is the total hamiltonian of the system. The standard form of a general dissipative 
process in Quantum Mechanics is governed by the Liouville equation obtained by adding a 
dissipation (super-)operator L D [p(t)} to Eq.®: 

ip(t) = [H,p(t)}+iL D [p(t)}. (9) 



|| For calculational simplicity we have chosen a degeneracy between the first and last energy levels, x\ = x±. 
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t {time in units of 1/y) 



Figure 1. Concurrence v Time t (units of 1 /control frequency) 



3.1. Liouville Dissipation 

In general, uncontrollable interactions of the system with its environment lead to two types of 
dissipation: phase decoherence (dephasing) and population relaxation. 

Phase decoherence occurs when the interaction with the environment destroys the phase 
correlations between states, which leads to changes in the off-diagonal elements of the density 
matrix: 

Pkn(t) = -i([H,p(t)))kn ~ ^knPknit) (10) 

where T kn (for k ^ n) is the dephasing rate between \k) and \n). 

Population relaxation occurs, for instance, when a quantum particle in state |n) 
spontaneously emits a photon and moves to another quantum state \k), which changes the 
populations according to 

Pnn 

[t) = -i([H,p(t)]) nn + ^2 [lnkPkk{t) ^JknPnn 

(<)] (ID 

kj^n 

where ^knPnn is the population loss for level |n) due to transitions |n) — > \k), and jnkPkk 
is the population gain caused by transitions \k) — > \n). The population relaxation rate jkn 
is determined by the lifetime of the state \n), and for multiple decay pathways, the relative 
probability for the transition \n) — > \k). 

Phase decoherence and population relaxation lead to a dissipation super- operator 
(represented by an iV 2 x iV 2 matrix) whose non-zero elements are 

(LD)[k;n],[k;n] = ~^kn k ^ 71 

(L D )[n;n],[h;k] = +lnk k^U (12) 

[LD)[n;n],[n;n] = —J2n^klkn 

where Tk n and "fkn are taken to be positive numbers, with Tk n symmetric in its indices. 
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In Eq. (fT2l we have introduced the convenient notation [m;n] — (m — l)N + n. The 
N 2 x A^ 2 matrix super-operator L D may be thought of as acting on the N 2 -vector r obtained 
from p by 

r[m;n] = Pmn- (13) 

The resulting vector equation is 

i- = Lr = {L H + L D )r (14) 

where Lh is the anti-hermitian matrix corresponding to the hamiltonian H. We obtain Lh 
explicitly by using the standard algebraic trick applied in evaluating Liouville equations (see, 
for example [5]). The correspondence between p and r as given in Eq. (fT3l) tells us, after some 
manipulation of indices, that 

p -> r => ApB ->■ A g) 5 r (15) 
using the direct (Kronecker) product of matrices. 

4. Physical Processes 

The quantum Liouville equation © is very formal; it covers both physical and non-physical 
processes and may tell us little about an actual physical dissipation process. For example, the 
values of the dissipation parameters Y kn and 7fc„ are not determined and a general choice will 
not lead to a physical process - that is, one under which the state p{t) remains a physical state 
- unless the parameters satisfy various constraints J6||. This (completely) positive evolution 
and the appropriate constraints emerge from physical stochastic dissipation equations such 
as those given by Lindblad and others in differential form 0, as well as in global form [[8]|. 
Nevertheless, the virtue of Eq.© is that essentially every dissipation process will have to 
satisfy it and so results derived from its use will have great generality. 

Since in our examples we wish to restrict ourselves to physical processes, we obtain our 
dissipation parameters Y and 7 by use of the Lindblad equation. 

4.1. Lindblad Equation 

Completely positive evolution of the system is guaranteed by the Lindblad form of the 
dissipation super-operator L D 

1 N 2 

LdW)} = -J2{[V s p(t),V s t] + [V s ,p(t)V}}} (16) 

A s=l 

where the matrices V s are arbitrary. The standard basis for N x N matrices is given by 

(Eij) mn = 5 im 6 jn (i, j,m,n — 1...N) (17) 
Relabelling, using the notation [m; n] — (m — 1)N + n, we choose 

V [i;j) = a [i;j] E [i;j]- ( 18 ) 
All the T's and 7's are determined by the (absolute values of) the A^ 2 (=16 here) parameters 
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4.2. Pure decoherence only 

So far we have discussed the most general case, when in principle all relaxation and 
decoherence parameters may be present in the dissipation matrix. However, experimentally, 
the relaxation time T\ for most systems is much longer than the dephasing time T 2 so that 
we may effectively neglect the relaxation rates 7. In the pure decoherence {dephasing) case, 
comparison of Eq.(fT6l) and Eq.(fT8l) with Eq. (fT2l) tells us that the 7 terms vanish if we choose 
a[i-,j] = for i ^ j. The decoherence parameters are then given by 

r« = \{H-a\ 2 + \ a m\ 2 ) (i,j = i...iv i^j)- (19) 

This leads to a mathematically very simple situation, as the dissipation matrix Ld is 
then diagonal. For the N = 4 system, this gives 6 pure dephasing parameters (Ty = and 
Ya = 0), determined by 4 constants, so there are two relations between the Ps - see Eq.(l2TI) 
below. 

Explicitly for the two-qubit, 4-level case, 

l do = diag{o, -r i2 , -ri 3 , -r M , -r 2 i, o, -r 23 , -r 24 , 

- r 3 i, -r 32 , 0, -r 34 , -r 4 i, -r 42 , -r 43 , 0} (20) 

with Tij = Tji. 

The constraints imposed by the physical process are 

Ti2 + r 3 4 = "14 + r 23 = r 13 + r 24 . (21) 



4.3. Concurrence 

In Sections [3] and |4] we are perforce dealing with general states, and so we must use the 
extended definition of concurrence for such (mixed) states [9J: 

Definition 7 (Concurrence: General definition) The concurrence C of a (mixed) two-qubit 
state p is given by 

C = max{Ai-A 2 -A 3 -A 4 ,0}, (22) 
where the quantities Aj are the square roots of the eigenvalues of the 4x4 matrix 

p{pi ® (? 2 )p*{o2 ® er 2 ) (23) 

in descending order, where cr 2 = ^ ^ 

This applies whether the density matrix p is either pure or mixed; for pure states it reduces to 
the form given in Definition [2l As already implied in the footnote given after Definition [2l the 
entanglement of formation is a monotonic function of the concurrence C, varying between a 
minimum of zero for C = 0, and a maximum of 1 for C = 1. 



4.4. Decoherence of Bell State 

We now give an example of a standard dephasing process acting on a maximally entangled 
state. 

Example 5 Consider the Bell state vg = l/v2[l, 0, 0, 1]. The Liouville vector r 
corresponding to this is 1/2[1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1]. The action of the 
dephasing operator L D0 ofEq.&M is given, as in Eq.KJM. by 

r = Lr= (L m )r (24) 
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concurrence 




Figure 2. Concurrence v Time t (units of 1/Decoherence T14) 



which may be immediately integrated to give 

r(t) = 1/2[1, 0, 0, e" ri4i , 0, 0, 0, 0, 0, 0, 0, 0, e~ Tl4t , 0, 0, 1] 
corresponding to the density matrix 

1 e- Tl4t 



(25) 



p(t) = 1/2 






-Tut o 






1 



(26) 



Note that this does not represent a pure state except att = 0. The concurrence as defined in 
Definition\7\evaluates to exp(— r 14 t). (See Figure 2.) 

The results of Example \5\ are essentially unchanged in the presence of an additional 
free hamiltonian, since this commutes with the dissipation super-operator L D0 , and indeed 
commutes with Ln in general ifTOl . and only introduces a phase factor. 



4.5. Stable Bell state 

In general, entanglement will decay under the type of dissipative processes noted here. 
However, as is clear from the last example, under special values of the decoherence 
parameters, entanglement will be preserved. In the case of Example [5] when T 14 = 0(= T 41 ) 
then the maximal entanglement of the state does not decay. Of course, in general we are not 
able to specify the values of the dephasing T's, but one may predict theoretically which types 
of state will remain invariant under the appropriate decoherence parameters. 
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